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Freehand three-dimensional (3D) ultrasoimd is a flexible imaging technique which allows a 3D data set to be 
constructed of sequential B-scans from a conventional ultrasound scanner. Since the data is acquired over 
several seconds, physiological motion generates spatial artifacts in visualisations of the data. Consequently, an 
electrocardiogram (ECG) signal is often used to gate the acquisition of B-scans to a single point in the cardiac 
cycle. We present a technique which can remove temporal artifacts by using properties of tiie grey-scale B-scan 
data, obviating the need for an external gating signal. B-scans are acquired throughout the cardiac cycle, and 
any phase can be selected for subsequent visualisation. This enables limited real-time 4D display of the data. 
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ABSTRACT 

Freehand three'<Iimensional (3D) ultrasound is a flexible imaging 
technique which allows a 3D data set to be constructed of sequen- 
tial B-scans from a conventional ultrasound scanner. Since the 
data is acquired over several seconds, physiolog:vcal motion gener- 
ates spatial artifacts in visualisations of the data. Consequently, an 
clectrocar^ogram (ECG) signal is often used to gate the acquisi- 
tion of B-scans to a single point in the cardiac cycle. We present a 
technique which can remove temporal artifacts by using properties 
of the grey-scale B-scan data, obviating the need for an external 
gating signal. B-scans are acquired throughout the cardiac cycle, 
and any phase can be selected for subsequent visualisation. This 
enables limited real-time 4D display of the data. 

1. INTRODUCnON 

Three-dimensional (3D) ultrasound is a developing imaging modal- 
ity which enables better visualisation of structure and measuie- 
ment of volume than conventional 2D ultrasound. There are two 
differing approaches: a dedicated 3D volume probe can directly 
record a 3D data set, or multiple images (B-scans) from a conven- 
tional 2D probe swept across the volume of interest can be refor- 
matted to generate 3D data. 

If the anatomy being scanned is dynamic, e.g. the heart, or pul- 
sating arteries, then the time taken to acquire the data set has a crit- 
ical bearing on the quality of the data. E)edicatcd volume probes 
would seem more appropriate for such applications, since they can 
acquire many volumes per second 11,2), fast enough to freeze car- 
diac (and respiratory) motion during acquisition. However, vol- 
ume probes are still in their infancy compared to traditional 2D 
probes, which provide considerably better resolution with an asso- 
ciated improvement in diagnostic potential. 3D probes also restrict 
the size of the acquired volume. 

The alternative, known as freehand 3D ultrasound [3J, has no 
such restriction: the probe can be nK>ved in whatever manner is ap- 
propriate for the anatomy being scanned. Since the data is acquired 
from 2D probes, the resolution (at least within the plane of each of 
the B-scans) is very good. However, motion of the anatomy results 
in spatial distortion within the 3D data set> since each 'slice' of the 
volume is acquired at a slightly different time. Such distortion af- 
fects the clarity of 3D data sets {45, and also the reliability with 
which volume measurements can be made from the data 15). 

Anatotmcal motion falls into three broad classes. The effect 
of patient movement can be limited by attaching a coordinate ref- 
erence to the patient, near to the area under investigation [6]. The 
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location of the 2D probe is then recorded with respect to this, rather 
than an external, static, reference. Respiratory motion can be sup- 
pressed by always scanning during a single breath hold, which lim- 
its the practical acquisition time to roughly 20 seconds. Cardiac 
activity is the other main source of anatomical motion, and since 
this is (generally) repetitive, ECG gating can be used to ensure 
that B-scans used in any single visualisation originate from the 
same phase of the cardiac cycle [7, 8, 9, 10, 111. This gating can 
be applied pre- or post-acquisition; in the latter case, the ECG sig- 
nal and B-scan data are recorded throughout the cycle, elfectively 
generating an (interleaved) 4D data set. 

ECG gating is not without its disadvantages, though. Clearly, 
it requires additional equipment connected to the patient and an ex- 
tra data stream input to the computer (in addition to the 2D probe 
locations and the B-scan images from die ultrasound machine). 
All three data streams must be temporally calibrated with respect 
to each other, since the phase of the ECG may not correspond to 
die pulse in the anatomy being scanned. In some cases, for in- 
stance for foetal heart scans, an ECG signal is not available [121. 
Hence, alternatives to an ECG have been investigated. For Dopplcr 
ultrasound data sets, it is possible to use the Doppler information 
to gate the data [13]. However, for grey-scale data sets, we would 
either need to pre-sct the cardiac frequency using a Doppler pre- 
scan [14], or use a separate source for Doppler gating and for data 
acquisition [15]. Neither solution is ideal — the first makes an un- 
realistic assumption of precisely uniform cardiac motion, and the 
second requires simultaneous use of two ultrasound machines. 

We present an algorithm which can gate freehand 3D ultra- 
sound data, post-acquisition, based on informadon in die data it- 
self. This enables dynamic 4D display of data affected by car- 
diac motion, with significantly reduced temporal jitter artifacts, 
and without the need for any additional equipment A different 
technique, based on Fourier transforms of secdons of the entire 
image sequence, has been applied to foetal heart scans [121. In 
contrast, our technique allows for variations in the cardiac cycle 
during the scan, and can also be applied to 3D Doppler data sets. 
We show results here for (adult) cardiac and arterial data. 

2. DESCRIPTION 

Any cardiac signal derived from grey-scale ultrasound data will 
inevitably be noisy. Ultrasound data is itself noisy, and only a 
small part of each B-scan may be pulsing: other parts may even 
be pulsing with a different phase. In addition, movement of the 
probe as the volume is scaimed will generate changes in the B- 
scans (since the anatomy is different in each scan) which are more 
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Fig. 1. Estimation of the phase of the input signal. The signal (top» 
dashed) is Gist pass-band filtered (top* solid), then the dominant 
frequency is found. The responses to single period sine and cosine 
filters at this frequency (middle) are then used to estimate phase at 
each sample point (bottom). 



significant than changes caused by cardiac motion. It is therefore 
vital that we extract as much information as possible out of the 
grey-scale cardiac signal; we will look at this first, before going on 
to see how the signal can be derived. 

2.1. Extracting cardiac phase fhmi a noisy signal 

Figure 1 (top, dashed) shows a typical grey-scale cardiac signal, 
plotted against the number of the B-scan from which it was de- 
nved; a pulse is only barely apparent from this signal. We ex- 
pect (adtUt) cardiac signals to vary between 30 and 150 beats per 
minute, allowing for missed beats and extra systoles. Approxi- 
mately 25 B-scans are acquired each second, which implies a nor- 
malised frequency range (where I represents half the sample rate) 
from 0.04 to 0.2. We can therefore pass>band Alter the signal to re- 
move any frequencies clearly outside this range. The filtered signal 
is shown in Fig. 1 (top, solid), and the filter's frequency response 
in Fig. 2. Note that this filter has a linear phase response in or- 
der not to distort the signal within die pass-band, and a low ripple 
(0.5dB) to minimise frequency weighting within the pass-band. 

The next step is to find the dominant frequency, which we 
can do simply by using sine transfonns of the entire filtered sig- 
nal. Providing the data set contains multiple cardiac cycles, these 
have very shaip responses: an example is shown in Fig. 2 (bot- 
tom, solid). We cannot simply use this frequency to gate the input 
data, however^ since we expect some variation in cardiac motion 
throughout die data set. This can be due to slight irregularities in 
the heart beat» or the result of dropped frames during acquisition' . 

Having found the dominant frequency, wc can estimate the 
phase at alt points in the signal by filtering it with single period 



'All aljorithins are applied across acquired framet laiher than lime. 
Slight unsteadiness in ihc acquisition rate therefore resulu in ^iparent 
slight fi^uency shifts in the ugnaL 
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Fig. 2. Filter responses. A linear-phase pass-band filter (top) re- 
moves physically implausible cardiac frequencies. Narrow-band 
filters (bonom, solid) are used to determine the dominant fre- 
quency, which defines the broader response sine and cosine filters 
(bottom, dotted and dashed). 



sine and cosine filters at this frequency. The inverse tangent of the 
output of each of these filters gives an estimate of the phase of the 
dominant signal. Essentially the same technique can be used for 
decoding phase modulated signals. 

The sine and cosine fitters have a broader response than the 
firequency selection filter, as can be seen in Fig. 2 (bottom, dot- 
ted and dashed), and hence can accommodate variation in the fire- 
quency of the cardiac waveform. The ou^ut of each filter is shown 
in Fig. 1 (middle), and the resulting phase estimate in Fig. I (bot- 
tom). In providing the centre frequency for the phase filters, the 
frequency selection filter can be seen as optimising the amplitude 
response of the phase filters, and hence the robustness of the phase 
estimate. In all other respects, the phase estimate is entirely local 
and independent of ^ estimate of the dominant frequency. 

2.2. Generating a signal containing cardiac information 

Using this technique, we can estimate cardiac phase from signals 
containing limited cardiac information. We could potentially use a 
dcformable mesh ( 161 or image correlation [ 17] to track movement 
in regions of each B-scan. However, these approaches arc often ex- 
pensive in processing time, and arc limited by the effect of genuine 
anatomical changes between B-scans, as mentioned earlier 

In fact, we can adopt a much simpler approach: counting pix- 
els in each B-scan with intensities falling within a user-specified 
threshold. So long as we constrain the count to within the region 
of die image affected by cardiac motion (also user-specified), we 
can generate sufficient information to give a reliable estimate of 
cardiac phase with which to gate the data. 

Figure 3(a) shows one B-scan in a data sequence of the com- 
mon carotid aiteiy, with the user-specified circular region. The 
region is generated by simply clicking (to define the centre) tiicn 
dragging (to change the radius) in the B-scan. Such a region is 
defined in several B-scans within the sequence, as the location of 
the carotid artery with respea to the B-scan*s origin changes. The 
region used in any otiier B-scan is simply die linear interpolation 
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(a) Region of interest 



(b) Selected threshold 



(a) Region of interest 



(b) Selected threshold 



Fig. 3. Signal extraction from grey-scale data: conunon carotid 
artery. Image pixels are counted over a user-specified circular re- 
gion of interest (a) and grey-scale threshold (b). 



(c) Un-gated reslice 



(d) Gated reslice 



of the centres and radii of the nearest user-defined regions. The 
user then specifies an intensity range, shown grey in Rg, 3(b). The 
threshold is chosen to encompass a portion of anatomy whose area 
changes (at least within the selected region) with cardiac motion. 

Having defined the regions and threshold, the number of pix- 
els within the threshold and contained in the region are counted in 
each B-scan, and this provides the cardiac signal for subsequent 
phase esdmation. Phase can be estimated in just a few seconds, so 
the user can interactively define the threshold and region to opti- 
mise the phase estimate. The signal processing in Sec. 2.1 works 
equally well for Doppler data, for which the total number of red 
(moving) pixels can be used to provide the cardiac signal. 

Figure 4 (top) shows the display of the cardiac signal: having 
estimated this, the xiser determines the centre and range of phases 
to use in subsequent data visualisation. B-scans within this range 
are represented by vertical lines drawn on top of the cardiac signal. 
VisuaUsatioQ of the data, for instance the reslice in Fig. 4 (bottom), 
reacts dynamically to the phase sliders, allowing real-time display 
of the pulsating data. 

3. RESULTS 

Figure 4 shows a visualisation of a 3D data set of a carotid artery. 
The threshold, region of interest and resuldng signals arc as shown 
in Figs. 3 and K The pulsating artifacts in the un-gated data are 
clearly removed in die gated reslice in Fig. 4(b). where only one 
B-scan has been used from each cardiac cycle. However, the left 
hand side of the reslice is now poorly sampled compared to the 
original data: the probe was being moved too fast at this point. 

For reslices gated at maximum contraction or relaxation of 
the carotid artery, it is possible to use a group of B-scans for dae 
visualisation, rather than just one per cycle [18], since here the 
cross-^tional area is approximately constant for a short period of 
time. Figure 4(c) shows an example where B-scans with a range 
of phases have been selecutd. This clearly improves the resolution 
of the reslice, while still minimising the motion artifacts. 

The technique can also be applied to ultrasound scans of the 
heart, as demonstrated in Figure 5. Cardiac motion artifacts, clearly 
evident in the un-gated data in Fig. 5(c), are no longer apparent in 
the gated data. Fig. 5(d). 

Some care is required in the interpretation of the phase esti- 
mate in both cases, since it is derived fi-om the varying anatomical 



Tig. 5. Grey-scale gating in echocardiographic data. The region of 
interest in (a) and threshold in (b) have been used to gate a cardiac 
data set. A reslice through one of die ventricles (c) is dramatically 
in^roved by gating (d). 



data radier than an external ECG at a fixed location. The phase is 
therefore locked to the deviation of the data across all the recorded 
B-scans, rather than the actual phase of the cardiac cycle. A phase 
of zero therefore implies maximal contraction everywhere in the 
data sequence, rather than. say. the point of systole. 

Since our aim is to improve visualisation rather than estimate 
the cardiac cycle per se, this is not a problem. Indeed, it could be 
argued that for arteries, visualisation at the maximum contraction 
or relaxation all the way along the anery is ideal. However, the 
gated visualisation does noi strictly represent a freeze-frame of the 
moving anatomy, though in practice it is likely to be very similar, 
dependent on the setting of ti)e region of interest and threshold. 

4, CONCLUSION 

We have presented a technique which allows cardiac motion ar- 
tifacts to be removed from freehand 3D ultrasound images post- 
acquisition, without any additional equipment. The method re- 
quires only a few minutes of user interaction after data acquisition, 
and can be applied to either grey-scale or Doppler ultrasound data. 
Initial tests on arterial and cardiac data sets indicate the general 
applicability of this approach. The phase estimate is also capable 
of tracking variations in the cardiac cycle. Dynamic (40) (^splay 
of the data is possible, but care is required in the interpretation of 
such displays. 
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(c) Reslice with phase group 
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